APS/123-QED 



O 
O 



< 



O 

oo 
o 



- 1—^ 

X 



Comparative Study of Different Discretizations of the 4 Model 

Ishani Roy 1,2 , Sergey V. Dmitriev 3 , Panayotis G. Kevrekidis 1 and Avadh Saxena 2 
1 Department of Mathematics and Statistics, University of Massachusetts, Amherst, MA 01003-4515, USA 
2 Center for Nonlinear Studies and Theoretical Division, 
Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
Institute of Industrial Science, the University of Tokyo, Komaba, Meguro-ku, Tokyo 153-8505, Japan 

(Dated: February 8, 2008) 



3 



We examine various recently proposed discretizations of the well-known </> 4 field theory. We 
compare and contrast the properties of their fundamental solutions including the nature of their 
kink-type solitary waves and the spectral properties of the linearization around such waves. We 
study these features as a function of the lattice spacing h, as one deviates from the continuum 
limit of h — > 0. We then proceed to a more "stringent" comparison of the models, by discussing 
the scattering properties of a kink-antikink pair for the different discretizations. These collisions 
are well-known to possess properties that quite sensitively depend on the initial speed even at the 
continuum limit. We examine how typical model behaviors are modified in the presence (and as 
a function) of discreteness and attempt to extract qualitative trends and issue pertinent warnings 
about some of the surprising resulting properties. 



PACS numbers: 05.45.-a, 05.45.Yv, 63.20.-e 



I. INTRODUCTION 

In the past few years, a variety of physical applications ranging from Bose-Einstein condensates in optical lattices 
, to arrays of waveguides in nonlinear optics and even to the dynamics of DNA Q have stimulated an enormous 
growth in the study of discrete models and the differential-difference equations that describe them. These models 
have a two-fold role and importance. On the one hand, they serve as discretizations of the corresponding continuum 
field theories; however, on the other hand, they may also be important physical models in their own right, e.g. in the 
■ context of crystal lattices. 

' An important class of models (and discretizations thereof) that is relevant to a wide variety of applications consists 
\Q , of the so-called Klein-Gordon type equations 3 . One of the particularly interesting equations within this family is 
the so-called 4 model featuring a wave equation with a cubic nonlinear (odd-power) polynomial added to it. This 
model has been physically argued as being of relevance in describing cosmic domain walls in cosmological settings , 
but also structural phase transitions, uniaxial ferroelectrics or even simple polymeric chains; see e.g. [7| and references 
therein. A particularly intriguing feature that was discovered early on in the continuum limit was the existence of a 
fractal structure in the collisions between the fundamental nonlinear wave structures (a kink and an anti-kink) in 
this model. This is a topic that was initiated by the numerical investigations of Refs. 0,0 (see also and is a 
topic that is still under active investigation (see e.g., the very recent mathematical analysis of the relevant mechanism 
provided in Ref. 

On the other hand, more recently, an issue that has concerned research work has been how to produce discretizations 
of such continuum models (such as the 4 model or its complex cousin, the nonlinear-Schrodinger equation |lfj|) that 
preserve some of the important properties of the corresponding continuum limit. One of the non-trivial aspects 
of this endeavor is the generation of discrete models in space that maintain some of the key invariances of their 
continuum siblings. For instance, in the uniform continuum medium, solutions can be shifted arbitrarily along a 
certain direction x by any xq (x is the spatial coordinate and xq — const.), due to the underlying translational 
invariance. However, discretizations generically will fail to maintain that feature and in the most straightforward 
versions thereof, equilibrium static solutions exist for a discrete rather than a continuum set of xq .4]. Some of 
these equilibrium solutions correspond to energy maxima and are unstable, while others, corresponding to energy 
minima, are stable. The difference between such maxima and minima of the energy is typically referred to as the 
Peierls-Nabarro barrier (PNb) . One of the topics of intense research efforts in the past few years has been to develop 
discretizations that do not present such energetic barriers; this is done in the hope that the latter class of models 
may provide more faithful representations of their continuum counterparts, regarding both symmetry properties and 
traveling solution features. 

The result of the above considerations in recent years has been the systematic construction of a large class of 
non-integrable discrete Klein-Gordon equations free of the Peierls-Nabarro barrier (PNb-free). Since the standard 
discretization of the continuum (j) 4 models conserves the energy but does have a PNb, Speight and co-workers (see 
also 0,E1) originally used a Bogomol'nyi argument [l4|, in order to eliminate that barrier. A later successful attempt 
(that produced multiple PNb-free models) was based on a different perspective, namely the one of associating the 



2 



PNb-free models to momentum-conserving discretizations 15] . Subsequently, yet another such model was recently 
proposed in [l^. Furthermore, these approaches were systematized and generalized through their formulation by 
means of a two-point discrete version of the first integral of the static continuum Klein-Gordon equation [rH IT^ | . We 
should note in passing that similar discretization efforts have recently been extended to the nonlinear Schrodinger 
equation [Tl El l2lj . 

One of the important questions that naturally emerges in the presence of this extensive recent literature is, indeed, 
how accurately do we expect these models to track the continuum limit behavior and how various properties are 
affected by the discreteness, as a function of its characteristic parameter (the lattice spacing h). Already, to some 
extent, there have been concerns regarding that question in that simulations of more sensitive phenomena such as kink 
collisions in the "Speight discretization" |Iull3| were only faithful to the continuum limit for fine lattices (/ik 0.1 or 
less) |22j . Here, we examine this question a bit more broadly and in more detail through numerical computations and 
analytical considerations comparing/contrasting five different discretizations of the continuum <f> field theory. Among 
them is the "standard" , classical 4> 4 discretization (Model 1) 0, two energy-conserving discretizations, namely the 
Speight-one (Model 2) [ll| and the one of Ref. also labeled CKMS hereafter (Model 3), and two momentum- 
conserving discretizations stemming from Ref. [ija], labeled Kl (Model 4) and K2 (Model 5), respectively. For each 
one of these models, we begin by examining the properties of the fundamental building block nonlinear wave solutions, 
namely the discrete kinks (and antikinks). We show how to obtain such solutions analytically or semi- analytically 
and subsequently examine the spectrum of small-amplitude excitations (linearization) around them, among other 
properties (such as stability) because this spectrum plays a nontrivial role in the outcome of wave interactions. 
Finally, we focus on the latter (i.e., on solitary wave collisions between kinks and antikinks) and attempt to extract 
salient features of such interactions as a function of the lattice spacing h and for a set of different speeds (and also 
different initial separations between the kinks). 

Our main findings can be summarized as follows: 

• The different models yield kink profiles that are different between them and from the continuum limit. These 
differences are strongest for the CKMS model and are found to lead to shrinking kinks for the energy conserving 
models 1-3, while they lead to expanding kinks in models 4-5. The deviation from the relevant continuum profile 
grows as h 2 . 

• The boosting of the kinks in order to induce their collision excites their internal modes. This plays a significant 
role in the collisions, since for different initial distances, the excitation of the internal mode will carry a different 
phase (at the moment of collision) and may accordingly lead to different collision outcomes. 

• The different models have different properties as regards the elasticity of their collisions. The most inelastic 
collisions occur in the standard discretization of model 1. Perhaps the next least elastic collisions occur in the 
CKMS Model 3, then Kl (Model 4), Speight (Model 2) and K2 (Model 5) in order of increasing elasticity. 

• The elasticity of collisions changes as a function of the lattice spacing. In fact, remarkably so, the collisions are 
more elastic for larger values of the lattice spacing. This is also demonstrated in the decreasing dependence of 
the critical velocity (beyond which the solitary waves separate after one collision) as a function of h. 

The presentation of our results will be structured as follows. In section [H] we will present the various models and, 
in section IIIII compare their kink solutions and spectral properties. In section IIVI we will examine the properties 
of the collisions of the different models focusing on a few typical speeds of the incoming waves for different initial 
distances and for different initial lattice spacings. Finally, in section [V] we will summarize our findings and present 
our conclusions, as well as some open questions for future study. 



Starting from the continuum limit of the model, we note that the one-dimensional cf) 4 field is described by the 
Lagrangian C — Ek ~ Ep with the kinetic and potential energy functionals defined, respectively, by 



II. <j) 4 FIELD AND ITS VARIOUS DISCRETIZATIONS 




(1) 




(2) 
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where </>(x, t) is the scalar field of interest and subscript indices mean partial derivatives with respect to the corre- 
sponding variable. The resulting Euler-Lagrange equation is obtained by demanding that be a local extremum of 
the action S = J Cdt, and it reads 

fat = 0^ + 20(1 - (j?) . (3) 
The following kink (antikink) solution to Eq. @, 

0(s,t) = ±tanh j=±= , (4) 

is one of the simplest examples of topological solitons. In Eq. Q, v is the kink velocity and xq is its arbitrary initial 
position (signalling the translational invariance of the continuum model discussed in the previous section). 
The first integral of the static version of Eq. © , 

U(x) = 4>l-{l-4> 2 ) 2 = 0, (5) 

plays an important role in our considerations. The integration constant was set to zero in Eq. (JSJ, which is sufficient 
for obtaining the kink solutions. The first integral can also be taken in a modified form, e.g., as 

u(x) = ±(j) x - 1 + 2 = . (6) 

We study various lattice dynamical equations obtained by discretizing Eq. J3} on the lattice x — nh, where 
n = 0,±1,±2..., and h is the lattice spacing. The general form of the lattice equations studied herein is 

4> n = A 2 </>„ + F(0 n _i, 4> n , 4> n +l) 

= D(<j) n -i,<f> n ,<j> n+ i), (7) 



where 



and, in the continuum limit (h — > 0), 



A 2 ^„ = (<t>n-l - 2<fin + <t>n+l), (8) 



F(0„_ 1 ,0„,0„ +1 )^20(l-0 2 ). (9) 



As mentioned previously, of particular interest will be the lattices whose static solutions, satisfying the three-point 
static problem corresponding to Eq. J7J, 

D(<j) n -i,(f> n ,<f> n +i) = 0, (10) 
can be found from a reduced two-point problem of the form 

I7(&,_i,&,) = 0, (11) 

or of the form 



u{<t> n - x , (f> n ) = 0. (12) 

Equations pit and i|12f) are the discretized first integrals (DFIs) obtained by discretizing Eq. JSJ and Eq. J^J, 
respectively jl7| (see also 0). 

Lattices whose static solutions can be found from the two-point DFI are called translationally invariant because 
the equilibrium solution can be obtained iteratively from the nonlinear algebraic equation starting from an arbitrary 
admissible value 4> n . In other words, translationally invariant lattices support a continuum rather than a discrete set 
of equilibrium solutions parametrized, e.g., by the position of the nth particle, <j> n . 

Let us now calculate the PN barrier in the models whose static solutions can be found from a two-point DFI, i.e., 
in the translationally invariant lattices. We consider a continuum set of equilibrium solutions parameterized by the 
position of the nth particle, (f) n , in the range <j) n £ [<p^' ,(p^'\, assuming that all values of <j> n within this range are 
admissible. 
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The work done by the inter-particle and the background forces (originating from the discretized background poten- 
tial) to move (quasi-statically) the nth particle from the configuration <jP^ to the configuration tfi^ is 

W n = D((f> n _x,(j) n ,(j> n+ i)A(f) n , (13) 
and the total work performed to "transform" the whole chain from tfi^ to ip^ is 

oo 

w= w "- ( 14 ) 

n— — oo 

However, in the translationally invariant models, the availability of a path of equilibrium configurations allowing to 
transit from ip^ to (j)^ leads to D = and thus, W n — for all n. This, in turn, results in W — 0. This result 
suggests that there is no energy cost to transform quasi-statically one equilibrium solution into another through a 
continuous set of equilibrium solutions. In other words, the height of the Peierls-Nabarro barrier calculated along 
this path is zero. For Hamiltonian lattices, the total work is path-independent (and equal to the potential energy 
difference between the final and initial state) and we can claim that, in such dynamical lattices, the PN potential is 
zero. For non-Hamiltonian lattices the work is path-dependent and we can only claim the absence of the PN barrier 
along the path considered above (which, however, is a natural one). While there are mathematical subtleties as 
regards whether this notion yields zero PNb more generally for translationally invariant lattices, this is the definition 
of PNb-free models that will be used herein. A more detailed examination of the, admittedly interesting, pertinent 
topics is outside the scope of the present study focusing on the comparison between different discrete 4 models and 
will be delegated to a future publication. 

In the following, we consider various discrete <f> models reported in the literature describing their kink solutions, 
their spectra of small-amplitude vibrations around vacuum solutions (cf> n = ±1), and also the spectra of lattices 
containing one static kink, revealing the kink's internal vibrational modes. Physical quantities conserved by the 
lattices are given, if they exist. These results will be quite relevant also in the discussion of kink-antikink collision 
outcomes. 



A. Classical discretization (model 1) 



The "standard" discretization of Eq. is Q 

4> n = A 2 0„ + 20„(1 - </> 2 ), 

and this is the only lattice in this study that possesses a Peierls-Nabarro barrier. 
Model 1 conserves the Hamiltonian (total energy) 



Hi 



E 



,(,2 



li 2 



(15) 



(16) 



Static kink solutions in model 1 exist only for those waves centered at a lattice site (unstable) or in the middle 
between two neighboring lattice sites (stable). Solutions can be found by various numerical techniques. As a first 
approximation, one can adopt Eq. to write the following approximate static kink solution 



cf> n = ± tanh[/i(n - x )]. (17) 

One can use this ansatz in a fixed point scheme (such as a Newton method) for xq = or x = 1/2 (mod 1), to 
identify the exact discrete static solutions <j> n . 

Subsequently, introducing the ansatz <p n (t) = 0° + e n (t) (where 0° is an equilibrium solution and e n (t) is a small 
perturbation), we linearize Eq. I|15l) with respect to e„ and obtain the following equation: 

S n = A 2 e„ + 2e n - 6(0°) 2 e„. (18) 

For the small-amplitude phonons, e„ = exp(ifcn + itot), with frequency u> and wave number k, Eq. I|18|l is reduced to 
the following dispersion relation: 

- 2 = ^sin 2 Q)-2 + 6(^) 2 . (19) 
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FIG. 1: Model 1, Eq. 1151 : Frequencies of the kink's internal modes for different magnitudes of the discreteness parameter h. 
Results for the on-site (circles) and inter-site (dots) kinks. Two solid lines show the borders of the spectrum of vacuum, Eq. 
12U> . The on-site kink is unstable because the spectrum contains imaginary frequencies, while the inter-site kink is stable. All 
three internal models are below the phonon band. 

From Eq. I)19|l . the spectrum of the vacuum solution, = ±1, is 
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uj 2 = 4 + sin 2 | '- 



(20) 



The spectrum of the lattice when linearizing around a static kink is shown in Fig. ^ 



B. Energy-conserving model 2 

Here we use the following DFI obtained from Eq. (|12|l . 



_ , 4>n-<t>n-l 1 . 4>l-l + + 4> 2 n , . 

u 2 = ± 1 + - =0. (21) 

The Hamiltonian, H = Ek + Ep, denned by Eq. Q and Eq. (0), can be discretized as follows, 

which gives the equations of motion of the energy-conserving model after Speight ^lj (see also [Tflj). 

d d 

4>n = -U 2 {(f) n -l,(l)n)^—U2((j>n-l,(l)n)-U2((f>n,4 l n+l)Tr—U2(4>n,^n+l) 
0(f>n 0<p n 

= fi + y) A 2 ^ n + 2cj) n ~ i [24> z n + {<j) n + ^„_i) 3 + (</»„ + 0„+i) 3 ]. (23) 
It is clear that the static solutions to Eq. (|23|l can be found from the two-point problem, Eq. I)21H . We have 



<M±1 



n 3 V3 / , 9 . 6 , 3 



T^ ±J 2-V^ ± ^ + 7? + 4 ' (24) 



where one can take either the upper or the lower signs. The kink solution can be obtained iteratively from Eq. I|24|l . 
starting from any \(f> n \ < 1. For the on-site and inter-site kinks one should take for the initial value <p n = and 
4> n =3/h— y/3 + 9/h 2 : respectively. 



6 



2 5 



CO 




FIG. 2: Model 2, Eq. 1231 . Same as in Fig. Notice the difference of the mode closest to the origin which in this case remains 
at to 2 = contrary to what is the case for model 1. There are three internal modes. However, note that a fourth mode appears 
for h > 0.85. There is a zero mode for all values of h indicating the absence of a PN barrier. 



The equation of motion, Eq. (|23[) . linearized in the vicinity of an equilibrium solution <^ yields 



h 2 



I h — ) A 2 e„ + 2e r . 



- [Gi^YSn + 3« + <t>l^Y{en + £ n -l) + 3(0* + €+lY^n + £n+l) 

The spectrum of the vacuum solution, <^ = ±1, is 



2 „ A~ h • 2 
U> 2 = 4 + 4 — r^— sir I - 



/l 2 



On the other hand, the spectrum of linearization around a kink is shown in Fig. [2] 



(25) 



(26) 



C. Energy-conserving model 3 

We take the DFI, corresponding to Eq. J^J), in the form 

U3 = Tpi&n - 0n-l) 2 - (1 — 4>n-l4>n) 2 = 0. 

The equations of motion of the model of CKMS 0] , 

V _ Uz{(j)n,<j)n+l) — Uz{(j>n-l,<j)n) 



"n+1 



^_ X )(1-^) 



can be obtained from the Hamiltonian 
where the potential V(4> n ) is given by 



/; 2 



1 - l - k \n 



V(cj> r , 



V{<t> n ) = -- 2 Ui, h2 



(27) 



(28) 



(29) 



(30) 
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FIG. 3: Model 3, Eq. 1281 : Same as in Fig. The vacuum solution <j> n — ±1 is unstable for h > 1. The upper edge of the 
phonon spectrum lies above the scale of the figure. There are three internal modes. A fourth mode appears for h > 0.6. Note 
that there is a zero mode for all values of h indicating the absence of a PN barrier. 



The exact static kink (antikink) solution is ^(J 

4> n = ±tanh[/3/i(n — Xq)], tanh(/3/i) = h, (31) 

where xq is the arbitrary position of the solution. 

Alternatively, the kink solution can be found from Eq. (|27|) . We come to the iterative formula, 

K-i ± h (32) 



l ± /i0„_i ' 



where one can choose either the upper or the lower signs and one can interchange <p n and 4> n -i- To obtain a kink 
centered on a lattice site, one should use as a starting point the value 4> n = 0, while for a kink centered in the middle 
between two neighboring sites, (j) n = 1/h — \J\jh 2 — 1. 
The linearized equation of motion reads 

l+(^ 2 -3)«) 2 + /i 2 «) 4 , s 



The spectrum of vacuum solutions (f)® = ±1 i 



is 



1-h 2 h 2 \2 

On the other hand, the spectrum of the CKMS lattice with a kink is shown in Fig. 01 



D. Momentum-conserving model 4 

Discretizing Eq. JSJ as follows, 

1 + h 2 

U A = —T2-(<t>n ~ 0n-l) 2 ~ (1 - <Pn-lK? = 0, (35) 

we come to the model reported in the work of 0] (motivated by its corresponding, so-called Ablowitz-Ladik sibling 
discretization for the nonlinear Schrodinger equation |23j) 

- U±(<j) n ,<j) n+ i) - Ui(^) n -i,(j) n ) 
= 2 2 

<Pn+l — <Pn-l 

= A 2 n + (0 n+1 +0 n _ 1 )(l- ( /) 2 ). (36) 
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FIG. 4: Model 4, Eq. 1361 . Same as in Fig. There are three internal modes including the up" = mode indicating the 
absence of a PN barrier. 



This non-Hamiltonian PNb-free model conserves the momentum [ia | which has the form: 

Pi = ^ <An(<?Wi - 0n-l). (37) 
n 

The exact static kink (antikink) solution is 

4> n = ±tanh[/3/i(n — xq)], tanh(/3/i) = = (38) 

where xo is the arbitrary position of the solution. 

Alternatively, the kink solution can be found iteratively from 



4> n -i ± h/Vi + h 2 ^ 



i±^ B _i/i/\/iW 



where one can choose either the upper or the lower signs and one can interchange <f> n and 4> n ~i- To obtain the on-site 
(inter-site) kink one should use as initial value <f> n = (<f> n = vH- h 2 /h — 1/h). 

The equation of motion, Eq. (|36|l . linearized in the vicinity of an equilibrium solution 0„ assumes the form 

e„ = (1 + h 2 ) A 2 e n + 2e n - W>° ) 2 (e„-i + e„+i) - 2$ (^°_ x + < +1 )e„. (40) 

The spectrum of vacuum, <^>° = ±1, coincides with that of model 1, Eq. I|2(J|) . 

However, the spectrum of the linearization around a kink is different as shown in Fig. 0] 



E. Momentum-conserving model 5 

Discretizing Eq. JSJ as 

u 5 = -^(4>n - n -i) 2 - 1 + €-i + 4>l\(€-i + <&) - \€-x4>l = o, (4i) 

we obtain another momentum-conserving model of the type of (see also 01) 

7 U 5 ((j) n ,4> n+ i) - U 5 (4> n -i,4> n ) 

V n = 2 2 

<Pn+l — (Pn-l 

= A 2 0„ + *S±L+*2zi( 4 - - 20 2 - 2 +1 ). (42) 
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FIG. 5: Model 5, Eq. 1421 . Same as in Fig. A unique feature of this discretization is the presence of the vibrational modes 
lying above the phonon spectrum. There are three internal modes (including the zero mode) and three additional modes for 
higher values of h indicating the presence of modes above the phonon spectrum. 



This non-Hamiltonian PNb-free model conserves the momentum of Eq. I|37fl • 

Static solutions in this model can be found iteratively by solving the quartic Eq. (|41fl . 

The equation of motion, Eq. (|42|l , linearized in the vicinity of an equilibrium solution 0° is 

£*n = ^2£n + En— 1 + e n+l 

- e - 1 + £ " +1 [^ n _ x f + 2(0°) 2 + (0° +1 ) 2 ] • (43) 

The spectrum of vacuum is the same as for model 2, Eq. (|26(l . 

Furthermore, the spectrum of the linearization around a kink is shown in Fig. 



III. COMPARISON OF KINK PROPERTIES 
A. Spectra of vacuum and kink's internal modes 

We have presented the spectra for the classical </> 4 model (model 1) and for the four models free of the Peierls- 
Nabarro barrier (models 2-5). All models share the same continuum limit; that is why, for small h, their properties 
are close and they only start to deviate from each other, as h increases. 

If we divide the models in groups by the quantities they conserve, then models 1-3 belong to the energy-conserving 
group while models 4 and 5 conserve the momentum of Eq. (|37[l . 

Models 3 and 4 have the static solutions derived in jig and |24J. Comparing the DFIs of these models, Eq. I|27|) 
and Eq. I|35|) . we can see that the solutions for model 4 can be obtained from those for model 3 by substituting 
h — > hj\J\ + h 2 . Exact static kink solutions are given for model 3 by Eq. (|31|1 or Eq. i|32|) and for model 4 by Eq. 
(JHEJ or Eq. ffity. 

Exact static kink solutions for model 2 can be found iteratively from Eq. (|24|l . while the ones of model 5 can be 
obtained by solving the quartic Eq. H41[) . For Model 1, the full 3-point problem of Eq. (|l(Jfl needs to be solved. 

Comparing the spectra of the vacuum (band edges of the spectra are shown by solid lines in Figs. H-EJ; we n °te 
that: 

• Model 4 has the same spectrum of vacuum as the classical model 1, and the width of this spectrum vanishes 
only when h — ► oo. The vacuum solution is always stable because uj 2 > for any h. 
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FIG. 6: Difference between the continuum and the discrete static kink profiles in the five models at h — 0.15. The continuum 
static kink is given by Eq. Q with v = 0. The kink in model 3 has the largest deviation from the continuum kink profile. 
Kinks in the energy-conserving discrete models (1 to 3) have a width smaller than the continuum kink while for the momentum- 
conserving models 4 and 5, the situation is reversed. 



• Models 2 and 5 have the same spectrum of the vacuum. The width of the spectrum vanishes at h = 1. Close to 
this value of h, the phonon spectrum is narrow and hence potential phonon radiation (of a kink-like structure 
due to resonance of internal mode harmonics with the phonon band) is minimized. The vacuum solution is 
always stable because uj 2 > for any h. 

• Model 3 has an /i-dependent cubic term; that is why the lower boundary of the spectrum is also /i-dependent, 
while in all other models it is constant (uj 2 = 4). In this model, the vacuum solution is stable only for < h < 1. 

Subsequently, examining the spectra of lattices containing a static kink, we note that (frequencies of kink's internal 
modes are shown in Figs. H-GOby circles and dots): 

1. Models 2-5 are PNb-free because they have a zero-frequency mode, which is the, so-called, translational (or the 
Goldstone) mode of the kink. For model 1, the corresponding mode has a non-zero frequency (in fact, it depends 
on h as oc exp(— n 2 /h); see e.g., [25j). signalling the presence of the Peierls-Nabarro barrier. Given the form of 
its /i-dependence, for small h (< 0.4), even for model 1 this mode has a nearly zero frequency (see Fig. QJ. This 
is the weakly perturbed translational mode of the continuum <j> equation. 

2. For small h (< 0.4) and even for moderate h (< 0.8), for all five models, apart from the translational mode 
we have two kink internal modes lying below the phonon spectrum, one of them very close to the edge of the 
phonon spectrum (u> 2 = 4) and another one in the vicinity of u> 2 ~ 3 (i.e., the corresponding continuum limit 
of this mode 26]). Additional internal modes may emerge for large h, but we will focus on smaller values of h 
(i.e., for h < 0.5) in the collision results that follow, hence we do not discuss these further here. 

3. Model 5, in contrast to all other models, is the model with internal modes lying not only below but also above 
the phonon band. Such modes, in contrast to the previously reported ones below the band, are of the short-wave 
(staggered) type and, for this reason, their excitation (or lack thereof) can be sensitive to the position of the 
collision point with respect to the lattice. 



B. Static kink profile and kink boosting 

It is also of interest to compare the static kink profiles in the five discrete models (to examine the relevant deviations 
between them and the continuum limit from which they are derived) . In Fig. we present the difference between 
the static kink profiles of the different models and that of the continuum static kink, Eq. J2J . The lattice spacing is 
h = 0.15. It is clear that the kink in the CKMS model 3 has the largest deviation from the continuum kink profile. 
Kinks in the energy-conserving discrete models (1 to 3) have widths which are smaller than that of the continuum 
kink while for the momentum-conserving models 4 and 5 the situation is reversed. Furthermore, in Fig. \7\ we show 
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FIG. 7: The maximal difference between the continuum and the discrete static kink profiles shown as a function of the spacing 
h for the five models. The amplitude of the deviation from the continuum kink profile increases with the discreteness parameter 
as h 2 . 



0.044 




FIG. 8: Kinetic energy of a lattice containing a kink boosted at t = with v = 0.25 with the help of Eq. IJ440 for the five 
discrete models. The lattice spacing is h — 0.3. 



how this difference between continuum and discrete kink is amplified as h is increased. This is carried out by showing 
the maximal difference between the two kinks as an /i-dependent diagnostic which reveals that the relevant difference 
grows as h 2 as h increases. 

In order to induce collisions, the kink needs to be set into motion in the discrete system. This can be achieved in 
a variety of ways. Here we have used the most standard one, namely Lorentz boosting the static kink to speed v, 
according to the continuum ansatz: 

(/>(x, t) = ± tanh v ; . (44) 

VT — v 2 

The results presented previously about the static kink also have a direct bearing on the kink boosting. In Fig. 
[S] we show the kinetic energy of a lattice containing a kink boosted at t = with velocity v = 0.25 through Eq. 
()44p. Lines of different thickness show the results for the five discrete models. The lattice spacing is h = 0.3 in this 
figure (even though similar, yet less pronounced results have been obtained for smaller h; again the relevant trend 
is quadratic in h). Ex(t) oscillates with the frequency close to the kink internal mode frequency of to 2 ~ 3. At 
t = we have Ek = 0.043 and for the energy-conserving models 1 to 3, Ej((t) is below this value, while for the 
momentum-conserving models 4 and 5, it is above this value. This is in sync with the static results where it was 
shown that models 4 and 5 have a correction to the continuum kink profile of opposite sign than the models 1 to 3 
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TABLE I: Domain Sizes and Initial Kink Separations (in units of lattice constant) Examined 



domain size 


80 


80 


160 


separation 


14 


28 


28 



(see also Fig. Model 3 shows the largest amplitude of kinetic energy oscillations, again in agreement with the 
static results. 

We have investigated another boosting method that uses the dynamical solution of the form </>„(£) = 0° + vte n , 
where </>° is the static kink solution, e n is the normalized translational kink's internal mode corresponding to the 
multiple eigenvalue uo 2 — 0, and v is the amplitude that plays the role of kink's velocity. We found this method to 
be very good (internal modes were not excited for any h) for small v, as it should be, because the accuracy of the 
linearized equations of motion increases as the eigenmode amplitude decreases. However, for velocities of order of 
v ~ 0.1 the accuracy of this method is insufficient because it does not take into account the Lorentz correction of 
kink's width. The ansatz Eq. (|44ll takes into account this correction, but it does not take into account the discreteness 
of media and, hence, naturally it is less accurate for large h. 

The best results were obtained for the use of Eq. (|44|l together with the addition of the kink's internal mode with 
the amplitude chosen to compensate the excitation of such mode. In the present study we did not use this more 
complicated/hne tuned method. 

As a result of these considerations, even for relatively small h, the kink boosted employing Eq. (|44|l carries an 
internal mode of non- vanishing amplitude. This internal mode often plays a nontrivial role in determining the outcome 
of the collision in what follows. 



IV. COLLISION RESULTS 

A. Numerical findings for different lattice spacings, initial speeds and kink-antikink separations 

We have carried out a comparative study of kink collisions under different discretizations. Our results have been 
obtained for different domain sizes (i.e., lattice sizes) and with different initial separations detailed in Table [I] As 
illustrated above, all of our models share the same continuum limit. We have compared the scattering properties for 
four different (dimensionless) velocities v = 0.21, 0.225, 0.24 and 0.255, respectively presented in Tables Hll VI We chose 
these velocities motivated by their (continuum limit) phenomenology in the detailed examination of ,6j. In each of the 
Tables lillVI the collision results are shown with an increment of 0.025 in the lattice spacing h for each of the different 
selected initial separations and domain sizes. Some of the standard collision outcomes are highlighted for v — 0.255 
in Figs. 191131 The most typical cases are those of figures 19*1101 the former shows a bion formation (i.e., the kink and 
the antikink merge, forming an oscillatory, so-called bion state, and never separate thereafter), a behavior typical for 
sufficiently small speeds. The latter illustrates what is characterized as a "one-bounce" separation, a behavior typical 
for sufficiently large initial speeds. However, the delicate structure of collisions for an intermediate range of speeds 
may lead to additional fine structure including multiple bounces before the eventual separation of the two kinks, as 
illustrated in Figs. 1111131 

The general trend as displayed in the table ITT1 is that for velocity v = 0.21, when the kinks collide, they form a 
bion state. In the bound state the kink and the antikink are trapped by their mutual attraction. In the table we 
see that, for small lattice spacings, the behavior of the different models is similar (as is expected, given the common 
continuum limit); on the other hand, the dynamics starts to diversify between discretizations, as the spacing is 
increased. Remarkably so, for larger values of h, we observe that the collisions are more elastic and, in fact, typically 
result in a single bounce for sufficiently large h. For velocity v = 0.225 (Table ITO)) the kinks are in a two-bounce 
window for small lattice spacing h but the change in outcome with increasing h is rather drastic (especially since 
the two-bounce is a rather fine-tuned collision outcome, where the internal modes control the resonant transfer of 
energy from and back to its original kinetic form 0,0, El)- We see a similar trend for the case with velocity v — 0.24, 
whereby the kinks form a bion state in the continuum limit but for increasing h we observe multiple bounces. For 
higher h, the kinks in all five models collide quasi-elastically i.e., with a single bounce. For a kink velocity of v = 0.255 
(which is close to the critical velocity, above which only single bounce phenomena occur in the continuum) we see an 
additional quite interesting feature. In this case, the outcomes for different models are so sensitive that they may not 
even converge for very small values of h. 

Our results overall indicate that the elasticity of the collisions depends strongly on the lattice spacing as well as 
on the details of the particular discretization. The collisions appear to be more elastic for larger values of h, a 
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TABLE II: Outcome of Kink-Antikink Collisions for v = 0.21. 



Results for velocity 0.21 


h 


Campbell et.al. 


Speight 


CKMS 


Kl 


K2 




80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


0.025 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


0.05 


biou 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


biou 


bion 


bion 


0.075 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


0.1 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


0.125 


bion 


bion 


bion 


4 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


2 


2 


2 


0.15 


bion 


bion 


bion 


4 


4 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


4 


bion 


bion 


0.175 


bion 


bion 


bion 


bion 


bion 


3 


bion 


bion 


bion 


2 


2 


2 


bion 


bion 


bion 


0.2 


bion 


bion 


bion 


bion 


bion 


2 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


2 


2 


0.225 


bion 


bion 


bion 


2 


bion 


bion 


bion 


3 


3 


bion 


bion 


bion 


1 


1 


1 


0.25 


3 


bion 


bion 


1 


1 


1 


bion 


bion 


bion 


2 


2 


2 


1 


1 


1 


0.275 


bion 


bion 


bion 


1 


1 


1 


2 


1 


1 


bion 


bion 


bion 


1 


1 


1 


0.3 


2 


bion 


bion 


1 


1 


1 


bion 


1 


1 


bion 


2 


2 


1 


1 


1 


0.325 


2 


2 


2 


1 


1 


1 


2 


1 


1 


2 


1 


1 


1 


1 


1 



TABLE III: Outcome of Kink-Antikink Collisions for v = 0.225. 



Results for velocity 0.225 


h 


Campbell et.al. 


Speight 


CKMS 


Kl 


K2 




80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


0.025 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


0.05 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


0.075 


2 


2 


2 


3 


bion 


bion 


2 


2 


2 


2 


2 


2 


bion 


bion 


bion 


0.1 


2 


2 


2 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


0.125 


2 


bion 


bion 


2 


2 


2 


4 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


0.15 


bion 


bion 


bion 


bion 


2 


2 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


2 


2 


0.175 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


2 


2 


bion 


bion 


bion 


0.2 


bion 


bion 


bion 


2 


bion 


1 


bion 


bion 


bion 


bion 


bion 


bion 


1 


1 


1 


0.225 


3 


bion 


bion 


1 


1 


1 


1 


bion 


bion 


bion 


2 


2 


1 


1 


1 


0.25 


2 


2 


2 


1 


1 


1 


1 


bion 


bion 


bion 


2 


2 


1 


1 


1 


0.275 


bion 


bion 


bion 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.3 


2 


bion 


bion 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.325 


bion 


bion 


bion 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 



feature which seems to be counter-intuitive given that discreteness in this type of models is perceived as a source 
of dissipation of kinetic energy j^]. On the other hand, discreteness leads to the excitation of additional internal 
modes (see Figs. [21 EI and in particular) and hence, potentially, to more exotic dynamical outcomes of the collisions. 
Furthermore, as h increases the width of the phonon band decreases, hence potentially limiting the range of resonant 
modes and therefore the amount of radiated energy (see also the relevant discussion below). This particular feature 
(apparent collision elasticity increase as a function of h) would be certainly wortwhile of a separate and detailed 
theoretical investigation. From the general trends of our results, we also observe that the most inelastic collisions 
occur for model 1, as might be expected by the presence of the PN barrier in that model. Finally, one more note 
of caution worth making here concerns the disparity between the different model results even for small h. It is clear 
that such phenomena as the outcome of collisions depend strongly and sensitively on a variety of factors (including 
e.g., the internal mode excitations, the location of collision, etc.) to an extent that one should not expect identical 
collision outcomes among these models even rather close to the continuum limit (which the models share) . This is also 
a result that partially defies the conventional wisdom that would suggest that different collisional outcomes across 
discretizations might be expected only when the length scale of discreteness (the lattice spacing) becomes comparable 
to the size of the kink. 

In what follows, we briefly analyze one of the sources of the above mentioned sensitivity of the collision outcome, 
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TABLE IV: Outcome of Kink-Antikink Collisions for v = 0.24. 



Results for velocity 0.24 


h 


Campbell et.al. 


Speight 


CKMS 


Kl 


K2 




80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


0.025 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


0.05 


biou 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


0.075 


bion 


bion 


bion 


2 


2 


2 


bion 


3 


3 


2 


bion 


bion 


bion 


bion 


bion 


0.1 


bion 


bion 


bion 


3 


2 


2 


bion 


2 


2 


bion 


2 


2 


bion 


2 


2 


0.125 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


3 


3 


3 


2 


bion 


bion 


0.15 


bion 


2 


2 


2 


1 


1 


bion 


bion 


bion 


bion 


3 


3 


1 


1 


1 


0.175 


bion 


bion 


bion 


1 


1 


1 


bion 


bion 


bion 


bion 


bion 


bion 


1 


1 


1 


0.2 


2 


2 


2 


1 


1 


1 


2 


3 


3 


1 


2 


2 


1 


1 


1 


0.225 


2 


bion 


bion 


1 


1 


1 


bion 


3 


3 


1 


1 


1 


1 


1 


1 


0.25 


bion 


2 


2 


1 


1 


1 


bion 


bion 


bion 


1 


1 


1 


1 


1 


1 


0.275 


bion 


2 


2 


1 


1 


1 


2 


1 


1 


1 


1 


1 


1 


1 


1 


0.3 


bion 


bion 


bion 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.325 


2 


2 


2 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 



TABLE V: Outcome of Kink-Antikink Collisions for v = 0.255. 



Results for velocity 0.255 



h 


Campbell et.al. 


Speight 


CKMS 


Kl 


K2 




80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


80/14 


80/28 


160/28 


0.0125 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


bion 


3 


bion 


bion 


bion 


bion 


bion 


0.025 


bion 


bion 


bion 


2 


bion 


2 


bion 


bion 


bion 


bion 


bion 


bion 


2 


2 


2 


0.05 


2 


bion 


bion 


bion 


2 


2 


bion 


bion 


bion 


2 


2 


2 


2 


4 


4 


0.075 


bion 


bion 


bion 


1 


bion 


2 


bion 


2 


2 


bion 


2 


1 


1 


1 


1 


0.1 


2 


bion 


bion 


1 


1 


1 


1 


bion 


bion 


1 


1 


1 


1 


1 


1 


0.125 


bion 


bion 


bion 


1 


1 


1 


1 


bion 


bion 


1 


1 


1 


1 


1 


1 


0.15 


1 


3 


3 


1 


1 


1 


1 


2 


2 


1 


1 


1 


1 


1 


1 


0.175 


1 


bion 


4 


1 


1 


1 


1 


2 


2 


1 


1 


1 


1 


1 


1 


0.2 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.225 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.25 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.275 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.3 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.325 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 



namely the original distance between the kink-antikink pair. We also quantify in a characteristic, in our view, way the 
increase in collision elasticity through the dependence of the critical speed separating bion formation from one-bounce 
collisions as a function of the lattice spacing h. 

B. Initial distance between colliding kinks 

In Fig. El we show the kink velocity after the collision as a function of initial collision velocity in model 3 at 
h = 0.15. Dots show the results for initial distance between kinks of 14, while open circles for the initial distance of 
14.8. 

One can notice a strong sensitivity of the collision outcome to the initial separation distance. This is because of 
the internal mode being excited when boosting the kinks. Changing the initial distance, we change the phase of the 
internal mode at the collision point, which critically, in turn, affects the result of the collision. 

For different models, the sensitivity of the collision outcome to the initial kink separation correlates with the 
amplitude of the kink's internal mode excited at boosting (see Fig. |SJ. Thus, the sensitivity is highest for model 3 
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FIG. 9: Bion formation for model 3 with h = 0.05, kink velocity = 0.255, domain size = 80, kink separation = 28. Note the 
radiation (phonons) emanating from the bion. 



and lowest for model 5. 

The sensitivity also decreases rapidly with decrease in h and the reason is, essentially, the same: the amplitude of 
the excited kink's internal mode decreases as h 2 . 



C. Threshold velocity v c as a function of h 

It is well known (see, e.g., 0) that there exists a threshold velocity v c such that collision of kinks with v > v c leads 
to separation after the first collision while for collisions with v < v c , the first collision cannot lead to separation. In 
the latter case, the reflection windows discussed in can be observed amidst regions of bion formation. 

In Fig. 1151 we show for the five models how v c changes with h. These results were obtained from computations 
similar to the ones presented in Fig. EI v c was estimated from the fit constant x (v 2 — v 2 ) 1 / 2 suggested in Q and 
thus the effect of the initial separation was averaged out to some extent. 

In all models v c decreases with increase in h implying that for larger h the collisions are more elastic. The standard 
discretization (model 1) shows the weakest dependence of v c on h, while in models 2 and 5 this dependence is strongest. 

Since the PN barrier is very small at h ~ 0.1, the observed effect can hardly be explained through the influence of 
the PN barrier. As one possible explanation of the dependence of v c on h, we discuss the burst of radiation emitted 
during collision. Corresponding numerical results are presented in Fig. 1161 for kinks colliding with v = 0.26 in the 
Speight lattice (model 2) with two lattice spacings, h = 0.1 and h — 0.15. This is the highest velocity we use in our 
simulations. We show the kinetic energy of radiation, (in order to exclude the kinetic energy of the moving kinks, an 
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FIG. 10: One bounce for model 2 with h = 0.1, kink velocity = 0.255, domain size = 80, kink separation = 14. 



area of width equal to 4 around each kink was not included in the computation of the kinetic energy) as a function 
of the time after collision. Dots show the results for h = 0.1 and open circles for h = 0.15. The amount of radiated 
energy grows with time due to the emission from the kink's internal modes excited at the collision. 

Extrapolation of the data presented in Fig. 1161 to t = suggests that, in the case of h = 0.1, the collision results 
in the burst of kinetic energy (in dimensionless units) of 8.7 x 10 -3 , while a smaller burst of radiation of 6.8 x 10~ 3 
takes place in the lattice with higher discreteness of h = 0.15. The fact that the burst of radiation is smaller in the 
lattice with higher discreteness can be related to the phonon spectrum width, which decreases with h as 1/h for small 
h, for all five models. The narrower the phonon band, the smaller the amount of energy that can be radiated and the 
more elastic the collision. 



V. CONCLUSIONS AND FUTURE CHALLENGES 



In the present work, we have analyzed the properties of a number of recently proposed discretizations of the 
continuum </> 4 field theory in the vicinity of (and further away from) the continuum limit. The relevant analysis 
consisted of the examination of the static properties of the models, concerning their fundamental nonlinear wave 
solutions, namely the kinks (and anti-kinks). For these types of solutions, we have examined how to obtain them, in 
what ways they differ from their continuum siblings, as well as the spectral properties of the linearization around such 
solutions. In particular, we have computed both the phonon (continuous) spectrum, as well as discussed the internal 
(or shape) modes present in the models. 

On the other hand, we have also examined dynamic properties of the kinks by studying their collision and com- 
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20 40 60 80 100 120 140 160 180 200 

FIG. 11: Two bounces for model 4 with h = 0.05, kink velocity = 0.255, domain size = 80, kink separation = 28. 



paring/contrasting their outcomes across the different discretizations. In that regard, we have observed a variety of 
interesting results. In particular, we have seen that the models only align with their continuum limit (especially as 
regards sensitive collision phenomenology) extremely close to the continuum limit [(i.e., for spacings of O(10 -2 ))]. In 
fact, in some cases (e.g. for the CKMS model and v = 0.255), the results are not independent of factors such as lattice 
spacing and kink-antikink separation even for the smallest h used herein (h = 0.0125). This is rather remarkable 
given that the scale of the kinks themselves is considerably wider, hence one would not expect this result on the basis 
of length-scale competition. However, we have argued that this should be attributed to the (initial-boost induced) 
excitation of the internal modes of the kink whose coupling to the continuous spectrum sensitively affects the collision 
outcome, as has been substantiated previously jBiHlSIIl- This should operate as a significant note of caution to 
researchers conducting numerical experiments with these models in an attempt to describe their continuum limits. 

Furthermore, we have seen that the elasticity of collisions varies not only from model to model, but even with 
increasing spacing of the lattice. In particular, we have illustrated through our numerical observations that the 
most inelastic collisions take place in the model that does have a Peierls-Nabarro barrier, while PNb-free models 
feature more elastic collisions. Moreover, there is a very interesting (and worthwhile to investigate further, possibly 
theoretically as well) dependence of the critical speed for single-bounce collisions on the spacing h. In particular, v c 
rapidly decreases as a function of h, rendering coarser collisions more elastic. This should also be attributed to the 
spectral properties of the models and the decreasing width of the phonon band for increasing h, which activates fewer 
couplings of internal mode frequency harmonics with the continuous spectrum and hence leads to weaker "dissipation" 
and consequently to more elastic collisions. 

While the static properties of the kinks have been obtained to a large extent explicitly from the underlying discretized 
first integral formalisms, kink collisions are naturally much harder to analyze theoretically for the presented discrete 
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FIG. 12: Three bounces for model 1 with h = 0.15, kink velocity = 0.255, domain size = 80, kink separation = 28. 



models. However, some of the relevant features such as the dependence of v c on h may be, to a certain degree, tractable 
(see e.g. |9| and references therein). Hence, it would be particularly interesting to seek a deeper understanding of the 
features numerically observed herein and how these can be associated with the nature of the underlying discretized 
nonlinearity. Such studies are currently in progress. 
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FIG. 14: Kink velocity after collision as a function of collision velocity in model 3 at h — 0.15. Dots show the results for initial 
distance between kinks of 14, while open circles for the initial distance of 14.8. 
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FIG. 15: Critical velocity v c as a function of h for the five models. In all models, v c decreases with increase in h implying that 
for larger h collisions are more elastic. The classical discretization (model 1) shows the weakest dependence of v c on h, while 
in the models 2 and 5 this dependence is strongest. 
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FIG. 16: Kinetic energy of radiation as a function of the time after a collision in model 2. Dots show the results for the lattice 
spacing h — 0.1 and open circles for h = 0.15. The collision velocity is v — 0.26. Radiated energy growth with time is due to 
the coupling of the kink's internal modes with the phonon band. Extrapolation of the data to t = suggests that in the case 
of h — 0.1 collision results in the burst of (dimensionless) kinetic energy of 8.7 x 10~ 3 , while a smaller burst of radiation of 
6.8 x 10~ 3 takes place in a lattice with higher discreteness of h = 0.15. 
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